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We compute all dynamical spin-spin correlation functions for the spin-1/2 XX Z anisotropic 
Heisenberg model in the gapless antiferromagnetic regime, using numerical sums of exact deter- 
minant representations for form factors of spin operators on the lattice. Contributions from in- 
termediate states containing many particles and string (bound) states are included. We present 
modified determinant representations for the form factors valid in the general case with string solu- 
tions to the Bethe equations. Our results are such that the available sum rules are saturated to high 
precision. We Fourier transform our results back to real space, allowing us in particular to make a 
comparison with known exact formulas for equal-time correlation functions for small separations in 
zero field, and with predictions for the zero-field asymptotics from conformal field theory. 

I. INTRODUCTION 

In recent years, enormous progress has been made in the quest to overcome an important and long-standing 
limitation of the Bethe Ansatai framework in the theory of integrable models: the inability to compute correlation 
functions. Most efforts have been focused on the Heisenberg spin chain2iMi&, although similiar sets of results could 
in principle be obtained for other integrable models formulated via the Algebraic Bethe Ansata&. In the particular 
case of the XXZ model, matrix elements of any local operator between two Bethe states can now be written as 
determinants of matrices, whose elements are known analytic functions of the rapidities of the eigenstates involved^. 
Used in conjunction with formulas for eigenstate norms£, this yields exact expressions for form factors on the lattice, 
thereby permitting in principle the computation of dynamical correlation functions. 

In this paper, we concentrate on the anisotropic Heisenberg spin-1/2 chain in a magnetic field, with Hamiltonian 



JY 



h =jy: 



(1) 



and periodic boundary conditions. The anisotropy parameter A can take on any real value without destroying the 
applicability of the Bethe Ansatz or of the present method. We will however consider here only the quantum critical 
regime —1 < A < 1. 

Our interest lies in space- and time-dependent spin-spin correlation functions. Continuing recent work by two of 
usm, we numerically obtain the dynamical spin-spin structure factor, defined as the Fourier transform of the connected 
spin-spin correlation function 

S aa (q,u;) =jjJ2 / *e*-*<5;(i)S?,(0)) c (2) 

jj'=i J -°° 

where a = z, ±. This quantity is of fundamental importance in many settings, and is used for example in the quanti- 
. , tative description of inelastic neutron scattering data for quasi-one-dimensional systemsSiifiiiAiiS to more theoretical 
questions relating to quantum entanglement^. In any case, the main interest is based on the fact that the Heisenberg 
spin chain is a strongly coupled quantum system with fractionalized (spinon-likc and higher) excitations, for which it 
■ - - i is possible to obtain nonperturbative information. 

The use of Bethe Ansatz-based methods to obtain quantitative results on correlation functions, in view of the long 
history of integrability, is rather recent. In the case of zero field, much progress has been possible for the model in 
the gapped phased, and for the isotropic limit, for which both two-spinon^^ and four-spinoni& contributions were 
obtained analytically. The nonzero field problem remains intractable with this approach. However, summations over 
intermediate states can be performed numerically once the analytical expressions for the form factors are known. 
Using this strategy, the spin-spin correlations (both longitudinal and transverse) for the isotropic Heisenberg model 
in a field were studied at fixed momentum q in |17|. and the transverse ones for XXZ at zero field and at q — n in 
|l8j. Two-particle contributions to the longitudinal structure factor for the XXZ chain in a field at q — n/2 were 
also studied in [n| ■ Multiparticle contributions and extension to the full Brillouin zone was done in [8j . In these, it 
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is observed that summing over a relatively small subset of intermediate states is sufficient to obtain good precision, 
as demonstrated by using sum rules. The present paper has three main further objectives: to generalize the form 
factor determinant formulae to intermediate states containing complex rapidities (in the form of strings) present in 
the Bethe Ansatz for the XX Z chain, to study the general field dependence of the correlation functions, and to obtain 
the real space- and time-dependent correlators. 

The plan of the paper is as follows. We begin in section [D] by recalling a few basic facts concerning the Bethe 
Ansatz for the XX Z model, and the lattice form factors for single spin operators. In section ITTT1 we discuss string 
solutions to the Bethe equations for excited states, and obtain modified determinant representations for the form 
factors of local spin operators involving such states. Section llVl contains the results of our numerical evaluations of 
the form factor sums for dynamical spin-spin correlation functions for two example choices of anisotropy. In section 
Ivl we use our results to compute the space- and time-dependent correlators. This allows us in particular to obtain 
equal-time correlators for any lattice separation. Comparison is made with known exact results in zero field for small 
lattice separation. We also compare our results to predictions from conformal field theory for the asymptotics in zero 
field. Our conclusions and outlook are collected in section IVTl 



II. SETUP 



In this section, we briefly review all elements necessary for the computation of the spin-spin correlation functions of 
the anisotropic Heisenberg model. The exact solution through the Bethe Ansatz for model JD is well-known (see [rlll2fH| 
and references therein). The reference state is taken to be the state with all spins up, |0) = ®fLi \ T)i- Since the total 
magnetization commutes with the Hamiltonian, the Hilbcrt space separates into subspaces of fixed magnetization, 
determined from the number of reversed spins M. We take the number of sites N to be even, and 2M < N, the other 
sector being accessible through a change in the reference state. 

Eigenstates in each subspace are completely characterized for 2M < N by a set of rapidities {Xj}, j — 1,...,M, 
solution to the Bethe equations 



sinh(Aj + i(/2) 
sinh(A J - iC/2) 



N 



M 

n 



sinh(Aj — Afe + i() 
sinh(Aj - Ak - i() 



j = l,...,M 



(3) 



where A = cos (. In view of the periodicity of the sinh function in the complex plane, we can restrict the possible 
values that the rapidities can take to the strip — tt/2 < ImA < 7r/2, or alternately define an extended zone scheme in 
which A and A + iirZ are identified. 

A more practical version of the Bethe equations is obtained by writing them in logarithmic form, 
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M 

-Y 



at an 



k=l 
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Here, Ij are distinct half-integers which can be viewed as quantum numbers: each choice of a set {Ij}, j — l,...,M 
(with Ij defined mod(iV)) uniquely specifies a set of rapidities, and therefore an eigenstate. The energy of a state is 
given as a function of the rapidities by 



M 



sin 2 C 



cosh 2 A, — cos£ 
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whereas the momentum has a simple representation in terms of the quantum numbers, 
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Ij mod 2n. 
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1, M, and all excited states are in principle obtained from the 



The ground state is given by Ij 
different choices of sets {Ij}. 

To study dynamics, some ingredients have to be added to the Bethe Ansatz: the matrix elements of spin op- 
erators between eigenstates (form factors). In terms of form factors for the Fourier-transformed spin operators 
Sq = -4= 5Zj=i e%q:> Sji t ne structure factor J2J can be written as a sum 



S aa (q,uj) 
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over the whole set of contributing intermediate eigenstates For the longitudinal structure factor, this is the set 
of all states with the same number of overturned spins as the ground state in the chosen magnetization subsector, 
excluding the ground state itself (connected correlator). For the transverse structure factor S it is the set of all 
states with one less overturned spin. oj a is the energy difference of state \a) with the ground state. Each term in (J7J 
can be obtained^ as a product of determinants of specific matrices, which are fully determined for given bra and ket 
eigenstates by a knowledge of the corresponding sets of rapidities. 

The form factor between two eigenstates (with M reversed spins, and rapidities {fJ.}, {A}) for the S z operator is 
given by£ 



\(W\S z q \{\})\ 2 = j5 q , q{x} _ VM 
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-i |det[H({ M },{A»-2P(M,{A»]r 
|det *({,;}) det*({A})| 
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where the matrices H and P are defined as 



H a6 (M,{A}) 
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Similarly, the form factor for the S operator between eigenstates with M and M — 1 reversed spins (respectively 
having sets of rapidities given by {/i}, {A}) is^ 



Yl™ |sinh(/i,- - iC/2)| 2 " , , , , 

l(MI^|{A})| 2 =N6 q , qw - qi J S m{\ ^ , 7 J n |sinh 2 ( Mj - Mfe )+sin 2 C| 

n j= i IsmhCAj- - *C/2)| i>fc=1 
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in which the 77 matrix is defined as 
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sinh 2 [i a + sin 2 C/2' 

The norm of the eigenstates appears in the above expressions in the form of the determinant of the Gaudin matriji^Si 



*ab({A}) = Sab 
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III. STRINGS 



Although the ground state of the XX Z chain is characterized by real rapidities, this is not true for all excited states. 
It has been known ever since the original work of Bethei that there exist sets of complex rapidities satisfying the 
Bethe equations. Typically, complex rapidities come in conjugate pairs or more elaborate structures involving higher 
numbers of elements. These complex rapidity groupings represent bound states of overturned spins, with spatial 
extent inversely related to the imaginary parts. At low densities (so at high magnetic fields), most solutions take 
the form of so-called strings, in which a number of rapidities share a real center but have regularly-spaced imaginary 
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parts. For the XX Z mod el, a classification of the possible types of string solutions has been conjectured in |22|, and 
further justified in |23l24j . This classification depends sensitively on the value of the anisotropy parameter A, and 
string state counting is traditionally used to argue for the completeness of the string state basis both for the isotropic 
magneliSiiS as for the anisotropic one2&. 

Our principal aim in this section is to adapt the determinant formulae for form factors and state norms to the 
case where one of the eigenstates involved contains string solutions to the Bethe equations. This is a necessary 
procedure, as both the Bethe equations and the determinant formulae for form factors suffer from degenerate limits 
in the presence of strings, which have to be analytically dealt with by hand. Let us however begin by reviewing the 
string structure, and setting out our conventions and notations. 

Following Takahashi and Suzuki^, we take C/ 71 " to be a real number between and 1, which is expressed as a 
continued fraction of real positive integers as 

' 1 ' ! -A > 1,14 > 2. (14) 

For large N, rapidities congregate to form strings centered either on the real line (for positive parity strings) or on 
the axis irr/2 (for negative parity strings), 

A2"° = X2' +^\{n 3 + l-2a) + l^{l-v 3 )+^8 n a l' a 1 a = l,..., nj (15) 

where the allowable lengths rij and parities Vj — ±1 are to be determined. In a string configuration, the parameters 
da 3 ' a are exponentially suppressed with system size. 

The classification of allowable string types in the thermodynamic limit proceeds according to the following algorithm. 
First, the positive integer series y~i, yo, yi, yi and mo, mi, mi are defined as 

i 

y_i = 0, y = 1, yi = i>i, and yi = Vi-s + ^iVi-i, i = 2,...,l, m Q = 0, rrij = } i/ k . (16) 

k=l 

Lengths and parities are then given by (our conventions here have the advantage of giving a proper ordering of string 
lengths, rij > n k , j > k) 



Vi-i 



+ (j-mi)Vi, ^ = (-l) L(n '- 1)C/7rJ , m,< 3 <m l+l . (17) 



The total number of possible strings is N s = mi + 1, and the index j runs over the set 1, N s . The real parameters 
Xa J represent the centers of strings with length rij and parity Vj, and are hereafter noted as X 3 a , a = 1, Mj, where 
Mj is the number of strings of length rij in the eigenstate under consideration. We therefore have the constraint 
T,f=i njMj = M. 

In a string configuration, many factors appearing in the Bethe equations become of the indeterminate form 5/5. 
Remultiplying Q for each member of a particular string gets rid of these factors, and allows one to rewrite the whole 
set of Bethe equations in terms of a reduced set involving only the string centers Aj,. Doing this, one finds (for N 
even) the reduced set of Bethe- Takahashi equations^ 

ef (Xi) = E jk (\i-\ k p), (18) 

(k,l3)^(j,a) 

where 

, % sinh(A + i?(l — v.j) + m,C/2) - , , , , o , . o , , , , , , 

gj - -^ sin h(A + 4(l-4)-4/2) ' ^W=VrCWi-.l +2 W^U- 3 (^- + ,W. (19) 

For the state classification and computation, it is preferable to work with the logarithmic form 

N„ M k 

N6j(Xi) -EE ®M - A £) = 2np a ( 2 °) 

fe=l /3=1 

where P a is integer if Mj is odd, and half-integer if Mj is even. The dispersion kernels and scattering phases appearing 
here are 

0j(X) = 2vj atan [(tanrijC/2) - ^ tanhA] , 
Bjfe(A) = (1 - 5 njnk )9\ nj ^ nk \{\) + 26»| nj _„ fc | +2 (A) + ... + 28 nj+nk - 2 (X) + nj+ „ k (\). (21) 
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The energy and momentum of a string are given by 

"sinh(A£ + i\ (1 - vj) + %n£/2) 



K = -J rWTj Ti Po(A J Q )=«ln 

Vj cosniAa — cosnf 



sml 



(22) 



ih(A J Q +if(l-^)-w^C/2) 
so the total momentum is again expressible in terms of the string quantum numbers as 

q = nT v=+ + ^EE 7 « m0d2 ^ ( 23 ) 

j=l a=l 

in which T v=+ is the total number of excitations with positive parity. 

It is known that the string hypothesis is not valid in absolute generality, meaning that not all eigenstates are 
described by a Bethe wavefunction with rapidities arranged in string patterns. For example, in the sector with two 
down spins of the isotropic antiferromagnet, it was demonstrated^ that, for large lattice size N, there exist 0(y/N) 
pairs of extra real solutions replacing the corresponding expected two-strings (the total number of two-strings scaling 
like N). A fraction 1/vJV of two-string states in the two down spin sector thus fail to be captured by the string 
hypothesis. For the anisotropic XX Z model, again in the two down spin sector, it was argued^ that the number of 
violations of the string hypothesis remains finite for any A < 1 as N goes to infinity. Non-string states for the XX Z 
model were also studied in [2^. In both cases, however, the total number of solutions to the Bethe equations given on 
the basis of the string hypothesis correctly reflects the dimensionality of the Hilbert space. If a string goes missing, 
it is replaced by a pair of real rapidities. 

The problem of counting the number of violations of the string hypothesis for larger numbers of overturned spins 
remains generally open. The one limit in which analytic progress can be made is at zero field, where M = N/2. 
There, it is possible to show^Mi that, in the thermodynamic limit, the lowest-lying excited states involving complex 
rapidities can only take the form of cither complex conjugate pairs or self-conjugate quartets. Strings of length higher 
than two therefore seem to be excluded. We therefore expect to see a breakdown of higher-than-two strings as the 
magnetization decreases from its saturation value (i.e. the limit in which strings exist) down to zero. 

Since we do only a partial trace over the set of all intermediate states, including only a minuscule proportion of the 
whole set, we do not concern ourselves here with non-string states. This will have as a consequence that our results will 
be "at their worst" at zero field, where we only sum up states including two-string excitations, but will be essentially 
exact at higher fields. For each string state, we check explicitly that the deviations 5 in l|15fl are exponentially small, 
and therefore only include these contributions if the string hypothesis is explicitly verified. In any case, we observe in 
general that strings of length higher than two give negligible contributions to the correlation functions. We therefore 
use the string hypothesis here as a general basis for our state classification scheme. The quality of our results will be 
confirmed later on with the use of sum rules. 

Our terminology for the state classification is as follows. First, we define a base as the set {M^}, k = 1, ... t N a of 
numbers of each string type present in an eigenstate. Many eigenstates share a same base, and these are obtained 
by choosing different quantum number configurations {ij}, j = 1, —,N a , a = 1, ...,Mj. The base itself defines the 
requirements on these quantum numbers. First, P a must be integer if Mj is odd, and half-integer if Mj is even. 
Second, they should be contained in the set {— I^,, — + 1) — j-^oo}' where 



1 Ns 

iL = Lo-|JWi(°°) - E( M * - <W i*(°°)IJj ( 24 ) 

fc=l 

in which the notation [a\ j means that we take the highest integer or half- integer less than or equal to the argument 
a, depending on the parity of Mj. The base choice completely specifies the right-hand side of this equation, and 
therefore the limits 1^ . 

( 2P + 1 \ 

The total number of configurations sharing a given base is then given by Z?{m} = II j I °m . ) ■ The total number 

of possible string states is then upon first examination given by summing this over all possible bases, D = J2{m} D{m}- 
This line of reasoning is used to argue for the completeness of the string states. However, we find that within a 
given base, there exist a number of inadmissible configurations. These correspond to configurations that are parity- 
symmetric, {— I 3 a } = {1^} and that have a vanishing quantum number for at least one higher string, 1^ = with even 
rij, or that have two vanishing quantum numbers I 3 al Ip, j ^ k, with rij, odd. In the first case, there are rapidities 
exponentially close to iQ/2 (which represents an ill-defined limit when obtaining the Bethe equations), and in the 
second at least two exponentially close rapidities, meaning a wavefunction improperly defined by the straightforward 
Bethe Ansatz. The correct eigenstates should thus presumably be obtained by a proper limiting procedure. In 
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view of the correctness of the string hypothesis state counting—, these must correspond to non-string states. We 
will readdress the question of state counting and non-string states in a separate publication. For our purposes here, 
however, completeness is an academic question, since we only need to do partial traces over intermediate states (albeit 
including strings). 

The presence of strings in the distribution of the rapidities {A} of an eigenstate renders the determinants used in 
(JHJ) , (| 1 1 1) degenerate, by which we mean that columns of the H, P, H~ and $ matrices become equal to leading order 
in the string deviations <5. The matrix determinants must therefore be expanded to higher order, and column/row 
manipulations used to extract the first nonvanishing contribution. Overall, the form factors have proper limits. 

It is straightforward to show that the norm of the eigenstate with strings is given by the determinant of a reduced 
Gaudin matrix, 



rij — 1 

idet$({A})i= n n n K^ i,o -^ ,B+i )i [idct^i+o^ 

j\nj>l a a=l 

of dimension equal to the number of strings in the state {A}. Its matrix elements arc defined as 
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In the case where the eigenstate described by the set {A} contains strings and the one described by {/z} does not, 
we can write explicit regular formulas for the elements of the matrices appearing in the longitudinal and transverse 
form factors. We find that the H and P matrix determinants for the longitudinal form factor can be reduced, in the 
presence of strings, to the determinants of reduced matrices up to corrections of order of the string deviations 6, 

| det[if -2P]| = | det[H^ - 2P«]| + 0(6). (27) 

The matrix components of and Pi r \ for matrix elements where Xb is part of an n-string, are given by 
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where 
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(29) 



in which we have extended the definitions of Aj, in 1)15(1 to include a = and a = n + 1 elements for notational 
convenience. In the last formula above, the product is defined as unity for n = 1. Note that these remain M- 
dimensional matrices (the dimension doesn't change here, in contrast to the Gaudin case). A similar procedure can 
be employed to treat the other form factor: the first M — 1 columns of the reduced transverse form factor matrix 
H~i r ^ are given by the same expressions as for (bearing in mind that the set {A} then has M — 1 elements), and 
the M-th column remains unchanged from equation l|12|) in the presence of strings in the set {A}. 

Therefore, if {A} contains strings and {/i} does not, we find that the string deviations 6 cancel out in the product 
for the longitudinal form factor, and that (with exponentially small corrections of order S) we can write 



KMI^|{A})| S 
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M 
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sinh(^j — *C/2) 
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sinh(Aj 
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(30) 
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The transverse form factor similarly reduces to 



n A Ii |sinh(u 5 -iC/2)\ 2 " 

|(M|5-|{A})| 2 = NS q , [sinCl ^E^ , ^ 7-^ II Nnh^-^ + zOl^x 

n i= i |smh(A J -»C/2)| j¥fc=1 

tt 1 mi |detH-«({^},{A})| 2 , 
x TT sinh (A 5 -Afc + jC ) r- 1 rr~, V w ,m (31) 



det*({/i})det*W({A}) 



For proper string states, all the above functions and determinants are nondegenerate. In the next section, we use 
these results to numerically compute the form factors, and sum them to obtain the dynamical correlation functions. 



IV. DYNAMICAL SPIN-SPIN CORRELATIONS 



The strategy to follow is now clear. We compute the S zz and S + structure factors by directly summing the terms 
on the right-hand side of equation (JJJ over a judiciously chosen subset of eigenstates. The momentum delta functions 



are broadened to width e ~ 1/N using 6 e (x) 



1— t 



x / e in order to obtain smooth curves. 



We scan through the eigenstates in the following order. First, we observe that the form factors of the spin operators 
between the ground state and an eigenstate {A} are extremely rapidly decreasing functions of the number of holes 
that need to be inserted in the configuration of the lowest-energy state (in the same base) in order to obtain the 
configuration {/} corresponding to {A}. We therefore scan through all bases and configurations for increasing number 
of holes, starting from one-hole states for S zz , and zero-hole states for S h . Although the number of possible 
configurations for fixed base and number of holes is a rapidly increasing function of the number of holes, we find that 
the total contributions for fixed bases also rapidly decrease for increasing hole numbers. We therefore limit ourselves 
to states with up to three holes, corresponding to up to six-particle excitations. 

We can quantify the quality of the present computational method by evaluating the sum rules for the longitudinal 
and transverse form factors. Namely, by integrating over momentum and frequency, we should saturate the values 



^Lly S "{q,w) = --(S') 2 = - 



2M N , 
1-1 ) 2 



(32) 



In Fig. ^ we plot the longitudinal structure factor as a function of momentum and frequency for anisotropy A = 0.25, 
for four values of the magnetization. Fig. contains the transverse structure factor for the same anisotropy and 
magnetizations, whereas Figs. |21 and 0] give plots for another value of anisotropy, A = 0.75. 

For all intermediate states involving strings, we explicitly check that the deviations from the string hypothesis are 
small. We find in general that states involving strings of length higher that two are admissible solutions to the Bcthc 
equations for high enough magnetizations. At zero field, only two-string states have exponentially small deviations S, 
and all higher-string states must be discarded, consistent with [30.31]. 

The relative contributions to the structure factors from different bases is very much dependent on the system size, 
the anisotropy, and the magnetization. In general, we find that two- and four-particle contributions are sufficient to 
saturate well over 90% of the sum rules in all cases, for system sizes up to N = 200. Interestingly, however, we find 
that string states also contribute noticeably in many cases. For example, in Fig. we plot the zero-field transverse 
structure factor contributions coming from intermediate states with one string of length two and up to three holes. 
Around six or seven percent of the weight is accounted for by these states, and similar or somewhat lower figures are 
found in other cases. Strings of length higher than two do not contribute significantly. For example, we find only 
around 5.7e-8 % of the sum rule from states with one string of length three, for the longitudinal structure factor for 
A = 0.25 at M = N/4 with N = 128. For A = 0.75, we find 6.3e-7 %. For the transverse correlators, the figures are 
2.3e-12 % and 3.1e-12 %. Even though these numbers would increase if we could go to larger system sizes, we do not 
expect them to ever become numerically significant. 

The imperfect saturation of the sum rules that we obtain in general can be ascribed either to higher states in the 
hierarchy which are not included in our partial summations, or states that are in principle included, but which are 
rejected in view of their deviations from the string hypothesis. As the proportion of excluded string states to allowable 
ones can be rather large (ranging anywhere from zero to fifty percent), we believe the latter explanation to be the 
correct one. We will attempt to include these non-string states in a future improvement of our method. 
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FIG. 1: Longitudinal structure factor as a function of momentum q and frequency to, for A = 0.25, and M = N/8, N/4, 37V/8, 
and N/2. Here, N = 200 and all contributions up to two holes are taken into account. The sum rule is thereby saturated to 
98.6%, 97.2%, 97.0% and 99.8%. 

V. REAL SPACE- AND TIME-DEPENDENT CORRELATORS 

From our results covering the whole Brillouin zone and frequency space, it is straightforward to obtain space-time 
dependent correlation functions by inverse Fourier transform: 

(s; +1 (*)s?(o)) c = ±J2\( GS \ s V a )\ 2e ~ iqaj ~ iWat - ( 34 ) 

a 

In the present section, we compare these results to known exact results for equal-time correlators at short distance, 
and to the large-distance asymptotic form obtained from conformal field theory. This comparison can only be made 
at zero field, where both sets of results are known exactly. 

A. Small distances and equal time 

It is known that correlation functions of spin operators on the lattice can be represented in general as multidimen- 
sional integrals2^*^2i2i^iS. For example, the spin-spin correlation function at distance j is exactly represented by 
a j-fold integral. These are unfortunately rather difficult to evaluate cither analytically or numerically, and exact 
answers are only known for small j or in the asymptotic limit j — ► oo. Recently, large simplifications were obtained for 
integral representations of the emptiness formation probability^S^L^, and for third-neighbour correlations^ in the 
case of the i sotropic H eisenberg chain in zero field. Results for the anisotropic chain in zero field were subsequently 
obtained in |40l4ll43| . yielding exact results for equal-time correlators of spin operators on up to four adjacent sites. 
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FIG. 2: Longitudinal structure factor as a function of momentum q and frequency u, for A = 0.75, and M — N/8, N/4, 3N/8, 
and N/2. Here, N = 200 and all contributions up to two holes are taken into account. The sum rule is thereby saturated to 
98.6%, 97.0%, 95.4% and 97.8%. 



Explicitly, these correlators read 41 

(S?S? +1 ) = F++ - -, (SjSf +1 ) = F+~, (S^S] +2 ) = - 2F++ + -, (SJ S+ +2 ) = 2F -t+ ■ 



(S*S 



2F1 



2F_ 



—+- 



3Fi 



2F 



+++- 



2F 



+ H 



j u j+3/ — **'++++ — \~ + 4' $/+2/ — -' - h ■ 



(35) 



where all blocks F, representing correlation functions of 2 x 2 fundamental matrices composed of vanishing elements 
except for one, are given explicitly as polynomials of the two one-dimensional integrals 



f"/ 2 . 1 coshr^ 

Qn\J) = I UX : , 

>-»ir/2 smn x sinh J -qx 



r-" /2 , 1 d coshryx 
-oo-in/2 smh x d-q sinh 3 -qx 



(36) 



The explicit form of these polynomials can be found for the XX Z case in the second appendix of paper |41|. We 
compare these results to our own values, obtained from Fourier transform of the form factors, in tables HI and ITU The 
deviations between the results are overall even smaller than can be expected on the basis of our slight deviations from 
perfect sum rule saturation. 

For the particular case of A = 0.5, exact results at zero field were also recently obtained up to distances of eight 
sites 4 ^. We compare these results to our form factor calculations in table Hill 



B. The scaling limit 



The continuum limit of the anisotropic XX Z chain in the gapless regime — 1 < A < 1 is described by a Gaus- 
sian conformal field theoryMi 4 ^. Expressions for the lattice Hamiltonian and for local spin operators as asymptotic 
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FIG. 3: Transverse structure factor as a function of momentum q and frequency u>, for A = 0.25, and M = N/8, N/4,SN/8, 
and N/2. Here, N = 200 and all contributions up to two holes are taken into account. The sum rule is thereby saturated to 
99.3%, 97.9%, 95.9% and 94.6%. 



TABLE I: Comparison of equal-time correlation functions {SjSj + i) at zero field for A = 0.25 for small distances I — 1,2,3. 
Subscript p refers to the exact polynomial representation, whereas // refers to our results obtained by summing form factors 
for all states up to three holes, thereby achieving saturation of the sum rule to 99.88 % (S zz ) and 95.61 % (S ). 



1 


QZZ 

Dp 


QZZ 


Sp 




1 


-0.113489 


-0.113337 


-0.316807 


-0.311455 


2 


0.0129789 


0.0129605 


0.180965 


0.180967 


3 


-0.0163964 


-0.0163965 


-0.152364 


-0.152466 



expansions of local scaling fields were obtained to subleading order in |4fl47j , whose results we use here as comparison. 
For large distance, the correlators explicitly read 



(Sj +l (t)S+(0)) c ~ (-1)' 



1 



B 



A/2 

(MT* { (M-)5- 

1 fi fl + l-\ U- 



A/2 



1 //+ L 



47rV+L [2 \l- l+J r/ + /_)f 



2 — ?7 (h 1 



D 



4(1-77) V/- 1+ 



(-l)%/4 



1 - 



(37) 



with l± = I ± t/e. To compare with our results, we use a conformal transformation to finite size, so that the chiral 
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FIG. 4: Transverse structure factor as a function of momentum g and frequency w, for A = 0.75, and M = iV/8, iV/4, 3iV/8, 
and iV/2. Here, iV = 200 and all contributions up to two holes are taken into account. The sum rule is thereby saturated to 
99.3%, 97.8%, 96.5% and 98.8%. 
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M/N .0.5 



: V 

0- 
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FIG. 5: The two-string contributions to the transverse structure factor at zero magnetic field, as a function of momentum q 
and frequency w, and for anisotropics A = 0.25 and 0.75. The density scale has been enhanced as compared to that used in 
the previous figures. Here, N = 200 and contributions up to three holes are taken into account. The sum rule contributions 
from these states is 7.2 % and 6.3 % respectively. 
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TABLE II: Comparison of equal-time correlation functions (SjSj + i) at zero field for A = 0.75 for small distances I = 1,2,3. 
Subscript p refers to the exact polynomial representation, whereas // refers to our results obtained by summing form factors 
for all states up to three holes, thereby achieving saturation of the sum rule to 98.73 % (S zz ) and 98.89 % (S~ + ). 



1 


QZZ 

Dp 


QZZ 

a ff 


S p ~*~ 


o-+ 
°.f.f 


1 


-0.136265 


-0.136093 


-0.305461 


-0.304906 


2 


0.043132 


0.043427 


0.140331 


0.140553 


3 


-0.035605 


-0.035911 


-0.117158 


-0.117201 



TABLE III: Comparison of equal-time longitudinal correlation function (SjS z +l ) at zero field for A = 0.5 for small distances 
I — 1, 8. Subscript e refers to the exact multiple integral representation of reference |4^|. whereas // refers to our results 
obtained by summing form factors for all states up to three holes, thereby achieving saturation of the sum rule to 99.2 %. 



1 


ST 


QZZ 

Z>.f.f 


1 


-0.125 


-0.123817 


2 


0.027344 


0.027141 


3 


-0.024475 


-0.024380 


4 


0.010994 


0.010955 


5 


-0.011084 


-0.011076 


6 


0.006248 


0.006246 


7 


-0.006567 


-0.006580 


8 


0.004153 


0.004163 



distances become l± — — sin(-^-(/ ± t)). Here, the anisotropy appears in the parameter rj = 1 — — , and the amplitudes 
A, A, A z , B, B, B z are known exactly for zero field, e.g. 

T( n ) 1V 



A 



2(1 -Vf 



2 v / ^r(^) J 



exp 



dt 



sinh(rjt) 



t \sinht cosh((l — rj)t) 



rje 



-21 



(38) 



(the explicit expressions for the other amplitudes, which we do not reproduce here, can be found in |47|). 

Tables IIVI and contain comparisons of the numerical values of the equal-time spin-spin correlation functions 
computed using the CFT aymptotic expansions and our form factor calculations. Fig. contains plots of the results 
as a function of lattice distance. Once again, the deviations are even smaller than should be expected, considering the 
imperfect sum rule saturation the the fact that these remain finite size comparisons. The only substantial deviations 
occur at small distances, where the conformal approach yields incorrect results (as a comparison to the data from the 
previous subsection shows). 



TABLE IV: Comparison of equal-time correlation functions (SjSj +l ) at zero field for A = 0.25 for large distances I. Subscript 
CFT refers to the scaling prediction, whereas // refers to our results obtained by summing form factors for iV = 200, for 
intermediate states including up to three holes, thereby achieving saturation of the sum rule to 99.88 % (S zz ) and 95.61 % 

(s-+). 



1 


QZZ 

JCFT 


QZZ 


o-+ 

°CFT 




10 


9.89624e-4 


9.89375e-4 


7.31376e-2 


7.39883e-2 


25 


-3.77195e-4 


-3.77096e-4 


-4.41637e-2 


-4.40758e-2 


40 


1.13604e-4 


1.13539e-4 


3.42939e-2 


3.43211e-2 


55 


-1.11476e-4 


-1.11424e-4 


-2.95934e-2 


-2.95775e-2 


70 


5.70852e-5 


5.70336e-5 


2.69545e-2 


2.69692e-2 


85 


-7.21563e-5 


-7.21058e-5 


-2.56504e-2 


-2.56405e-2 


100 


4.71343e-5 


4.70918e-5 


2.52109e-2 


2.52193e-2 
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TABLE V: Comparison of equal-time correlation functions (SjSj + i) at zero field for A = 0.75 for large distances I. Subscript 
CFT refers to the scaling prediction, whereas // refers to our results obtained by summing form factors for N = 200, for 
intermediate states including up to three holes, thereby achieving saturation of the sum rule to 98.73 % (S zz ) and 98.89 % 

(s-+). 



1 


QZZ 

a CFT 


QZZ 


Q- + 
°CFT 


c-+ 

a ff 


10 


7.10534e-3 


7.26298e-3 


4.22187e-2 


4.36235e-2 


25 


-2.49782e-3 


-2.51931e-3 


-2.24156e-2 


-2.21879e-2 


40 


1.39631e-3 


1.39725e-3 


1.57658e-2 


1.58537e-2 


55 


-1.03764e-3 


-1.03205e-3 


-1.31010e-2 


-1.30427e-2 


70 


8.26181e-4 


8.18235e-4 


1.14797e-2 


1.15143e-2 


85 


-7.56746e-4 


-7.47886e-4 


-1.08240e-2 


-1.07873e-2 


100 


7.13747e-4 


7.04615e-4 


1.05099e-2 


1.05363e-2 




FIG. 6: Equal-time spin-spin correlation functions as a function of lattice distance, for anisotropies A = 0.25 and A = 0.75, 
for chains of N = 200 sites. Our form factor results (for intermediate states containing up to three holes) are compared to the 
asymptotics from CFT adapted to finite size. The deviations are well within those expected from the sum rules and finite size 
effects for the form factors, and subleading corrections to the CFT result. Example numbers are given in tables HVI and IVI 

VI. CONCLUSION 

The existence of determinant representations for lattice form factors of local spin operators makes it possible to 
compute dynamical correlation functions to a high degree of accuracy, and we have here attempted to push the results 
of |8| further by including the remaining spin-spin correlators, studying their general field dependence, including 
contributions from intermediate states containing complex solutions to the Bethe equations in the presence of strings, 
and computing explicitly real space- and time-dependent correlators. Although intermediate states with strings of 
length two do indeed give non-negligible contributions to the sum rules, we have found that strings of any higher 
length have form factors which are much too small to be of significance. 

We have demonstrated the reliability of our results by performing comparisons with exact integral representations for 
equal time, small-distance correlators at zero field, and exact formulae for the large-distance asymptotics (again at zero 
field). We do wish to emphasize that the present method, unlike these calculations, is not limited to zero field. Also, 
unlike other numerical methods based on the Density Matrix Renormalization Group (DMRG), it is by construction 
directly applicable to dynamics. In other words, we are able to numerically compute all spin-spin correlators either 
as a function of momentum and frequency, or space and time separation, for any value of magnetization, for system 
sizes well beyond the reach of any other numerical method. 
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